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ABSTRACT 

We present new Green Bank Telescope (GBT) 21 cm neutral hydrogen (H I) 
observations of a complete distance limited sample of 22 Hickson Compact 
Groups (HCGs) with at least four true members. We detected an average HI 
mass of 8 x 10 9 M (median = 6 x 10 9 M ), which is significantly larger than 
previous single-dish measurements. Consequently, the Hi-deficiencies for these 
HCGs have been reduced, although not completely eliminated. Spectral compar- 
ison of the GBT data with complementary Very Large Array (VLA) data shows 
significant H I excess in the GBT spectra. The observed excess is primarily due to 
the high surface brightness sensitivity of the GBT detecting diffuse, low column 
density HI in these groups. The excess gas forms a faint diffused neutral medium 
which is an intermediate stage in the evolution of high-surface brightness H I tidal 
debris in the intra-group medium (IGM) before it is fully ionized. The excess 
gas mass fraction, (M(HI)gbt - M(HI)vla)/M(HI)g.bt, for our complete sample 
varies from 5% to 81% with an average of 36% (median = 30%). The excess gas 
mass fraction is highest in slightly H I deficient groups where the tidal debris has 
had enough time to evolve. We also find the e xcess gas content increases with th e 



evolutionary phase of the group described in IVerdes— Montenegro et al.l (120011 ). 
Theoretical calculations indicate that an H I cloud of radius > 200 pc would 
survive in an IGM of 2 x 10 6 K for more than the typical dynamical lifetime of 
a group. However, smaller clouds get evaporated and assimilated into the hot 
IGM in a much shorter timescale. 



Subject headings: galaxies: evolution — galaxies: interactions 



-2- 



INTRODUCTION 



Mos t galaxies in the universe are not isolated; instead they are in groups or clusters of 
galaxies ( lTullylll987l ). Galaxy groups are an environment where galaxies evolve continuously 
as they interact with each other and the group potential. Grou ps also contain a significant 
fraction of their baryonic conten t outside the individual galaxies ( IFukugita. Hogan &: Peebles 
19981 ; IFukugita fc Peebles! I2006T) . The detection of a substantial intragroup medium (IGM) 
in X-ray (see lMulchaevll2000h provides a conclusive evidence of baryons associated with the 
group potential. There is also growing evidence that significant changes in galaxy morphol- 
ogy and truncation of star formation take place in groups and cluster outskirts, rather than in 



dense cluster cores (lDresslerlll980l ; iLewis et al.ll2002t iHelsdon fc Ponmanl 120031 ; ICapak et al. 



20071 . and references therein). Therefore, a full accounting of the gaseous components in 
galaxy groups and their influence on the evolution of member galaxies should offer an im- 
portant insight into the evolution of galaxies in general. 

Hickson Compact Groups (HCGs) are unique laboratories for studying the effects of 
multi-galaxy tidal interactions on the morphological and dynamical evolution of galaxies as 
well as the IGM. HCGs are groups of galaxies sho wing extreme p rojected densities and were 
identified using the following set selection criteria (jHicksonlll982l ): (1) four or more galaxies 
within 3 magnitudes of the brightest; (2) no other galaxy in the magnitude range of the 
group members between the group radiu^] Rq and 3i?c; and (3) mean surface brightness 
of the group, p,c < 26.0 mag arcsec -2 in the red (E) plate of the Palomar Observatory Sky 
Survey (POSS). A total of 100 compact groups w i th pr ojected galaxy density ranging from 
300 to 10 8 h~ 2 Mpc~ 2 were cataloged by [Hickson] ( 119821 ). The projected galaxy densities in 
these groups are among the highest in the local universe . Despite the theoretic al arguments 
that these groups cannot last long enough to be real ( IHernquist et al.lll995[ ). subsequent 
spectroscopic surveys have s hown that the majo rity are physical groups with low velocity 
dispersions of ~250 km s _1 ( iHickson et al.lll992l ). 



Morphological analysis has shown that compact groups contain a significantly smaller 
fraction of late-type (spiral and irr egular) galaxies as compared to that seen in the field 



f lRood fe William 



989c lKindllll990l ). Moreover, spiral galaxies associated with these groups 



tend to be deficient § in neutral atomic hydrogen (H I), containing on an a yerage only 40% of 



the expected H I based on their optical morphology and luminosity types ( 


Williams & Rood 




1987; 


Williams & Van Gorkom 


1995; 


Huchtmeier 


1997; 


Verdes— Montenegro et alJl2001 


^Stevens et al. 



1 Rg is the group radius defined by the smallest circle containing the centers of the group members. 

2 This is referred to as "H I-deficiency" and is defined as Defui = log[M(HI) pre d] — log[M(HI) 6 S ] (see 
Havnes fe Giovanellilll984 ). 
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2004 . and references th erein) . These galaxies are also de ficient in CO emission indicating low 
molecular gas content ( IVerdes— Montenegro et al.lll998l ). Tidal gas stripping, shock heating, 
and gas exhaustion by tidally induced star formation are some of the possible explana- 
tions for the observed def iciencies. H I deficient groups te nd to be brighter in X-rays than 



the non-deficient groups (IVerdes— Montenegro et al.ll200ll . VM2001 hereafter). This sug 



gests that the hot IGM facilitates the ionization of tidally stripped cold interstellar medium 
(ISM). However, the IGM is not always detected in the X-ray, especially in spiral dominated 
groups, possibly because of a smaller halo mass thus resulting in a lower vi r ial tem perature 
( IMulchaey et al.l Il996at iRasmussen et al.ll2008l ). For instance, lAracil et al.l ( 120061 ) detected 
cooler neutral gas (< 10 5 K) in the form of a cluster of 13 H I Lyman a absorption lines spread 
over 1000 km s _1 in the IGM of a galaxy groups at z = 0.0635 towards quasar HS 0624+6907. 



Atomic gas studies using the Very Large Array (VLA) of HCGs f 


Verdes—! 


VIontenegro et al. 


2001; 


Williams k. Van Gor 


com 


1995; 


Verdes— Montenegro et al. 


2005 


2002; 


del Olmo et al. 




2003: 


Durbala et al. 


2008: 


Williams et al. 


2002 


) provide insights on the fate of cold gas in 





HCGs, including ubiquitous tidal features and off-disk gas. Tidal interactions are important 
in all compact groups, and they produce complex, continuously changing gaseous and stellar 
structures. Frequent, multiple tidal encounters within the group potential may rapidly dis- 
perse H I tidal debris to a diffuse, extended medium. Cold ISM content dictates the current 
and future star forma tion activity in galaxies. Examining the Spitzer infrared imaging data, 
Johnson et al.l (120071 ) found that the most gas-rich groups also tend to be the most actively 
star forming. They also found evidence of rapid evolution of galaxy properties in response 
to the dynamical effects and the level of these activities directly impacts the morphological 
and spectral evolution of the galaxies. Diffuse emission features with broad (> 600 kms -1 ) 
line-widths may escape detection when observed with a relatively limited bandwidth and 
surface brightness sensitivity. The main aim of this work is to explore whether there exists 
significant diffuse H I emission previously undetected by the VLA and to determine their 
quantity, location, and nature. 

We have obtained high quality H I spectra for a complete distance limited sample of 
22 HCGs along with 4 additional sources using National Radio Astronomy Obsrevatory's 
(NRAofl Robert C. Byrd Green Bank Telescope (GBT). The excellent sensitivity and spec- 
tral resolution of the GBT provides the opportunity to study the physical state and the 
fate of missing cold gas in HCGs. The work presented here is part of a broader study of 
evolution of galaxies in compact group environment through the study of morphology and 
tidal signatures traced in the 21cm H I line. 



3 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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A description of our sample along with the details of our GBT observations are presented 
in § 2 and the results are presented in § 3. We discuss our findings and use comparative 
studies of the new GBT measurements and the previous VLA measurements to investigate 
the nature and the distribution of the neutral gas in HCGs in § 4. Finally, a conclusion 
is presented in § 5, summarizing the results and inferences drawn from this work. For this 
study, we have adopted the value of the Hubble constant, H = 75 km s -1 Mpc -1 . 



2. OBSERVATIONS 



2.1. Sample 



Identified by visual inspection of the Palomar Sky Survey plates and without any redshift 
information, the original Hickson catalog is fraught with projection effects and selection 
biases. Further examination of the catalog incorporating the redshift information has shown 
that about 1/3 of the gr oups are false groups, either in pro jection or part of a more massive 
group or a cluster (see iHickson et al.lll992l ; ISulentid 119971 ) . Understanding that HCGs are 
potentially diverse and heterogeneous entities, a new coherent sample of true groups from 
the Hickson catalog has been constructed for a meaningful statistical analysis. The criteria 
we used to define the new complete sample are: 



• having four or more true member galaxies meeting the original Hickson definition, 
excluding false groups by utilizing redshift information; 

• containing at least one spiral galaxy, to determine H I deficiency meaningfully; and 

• located within a distance of 100 megaparsec (Mpc) so that they can be studied with 
excellent sensitivity. 

These criteria yielded 22 groups which we will refer hereafter as the complete sample: 
HCG 7, 10, 15, 16, 23, 25, 30, 31, 37, 40, 44, 58, 67, 68, 79, 88, 90, 91, 92, 93, 97, and 100. 
We have also included four more groups that do not fit the sample criterion but had previous 
single dish and/or VLA data - HCG 18, 26, 35, and 48 (referred to as the "additional sources" 
hereafter). HCG 18 and 48 have less than four group members while HCG 26 and 35 are 
further than 100 Mpc. Therefore, these four groups will not be included in our statistical 
studies presented in § 14.41 The redshift, size and morphological information of each of these 
groups and their members are provided in Table [TJ 



Whether HCGs are r eal groups or chance pr ojections along radially oriented filaments in 
the large-scale structure (IHernquist et al.lll995l ) is a highly debated topic. An independent 
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identification of galaxy groups was performed by lYang et al.l (120071 ) using the halo occupancy 
model on the galaxy catalogs from the Sloan Digital Sky Survey (SDSS) and the Two Degree 
Field (2dF) galaxy redshift survey. Six of our groups, HCG 7, 16, 25, 35, 55, and 88, fall 
within the sky and r edshift cover age of these surveys and all of them have been identified 
as galaxy g r oups by lYang et al.l The SDSS a nd 2dF surveys are deeper than that used 
by iHicksonl (119821 ) and hence the lYang et al.l catalog includes additional fainter galaxies 
associated with these groups. They derived halo masses for these groups, which range 
between (2 - 115) x 10 12 M . The independent identification of all 6 HCGs within the SDSS 
and 2dF sky coverage lends further support that our sample refinement described above 
successfully identified real physical groups. 



2.2. GBT H I and OH Observations 

We conducted GBT observations of the 21 cm H I and 18 cm OH transitions for 26 
HCGs in May and October 2005 as part of the GBT05A-051 and GBT05C-005 programs. 
The data on HCG 67 and HCG 79 were severely corrupted by solar interference, and new 
observations were obtained in August 2006 as part of the GBT06B-053 program. We used 
the dual polarization L-band system that covers the frequency range between 1.15 and 1.73 
Gigahertz (GHz). Using 9-level sampling and two IF settings, each with 12.5 MHz total 
bandwidth, the spectrometer gives 8192 channels at 1.5 kHz (0.3 km s _1 ) resolution, covering 
a total velocity range of 2500 km s _1 . We employed the standard position-switching scheme 
by cycling through the ON-OFF sequence, dwelling for 300 seconds at each position. A 
position offset of +20' in Right Ascension was adopted for the OFF position so that the 
presence of any confusing sources in the OFF position could be tracked. We also performed 
local pointing corrections (LPC) using the observing procedure AutoPeak, which was then 
automatically applied to the data. Hence, the error in pointing should be less than GBT's 
global pointing accuracy of 2.7 ", which is less that 0.5% of GBT beam size at 1.4 GHz. 

The clean optics and active surface system of the GBT also results in a well-calibrated 
structure and a stable gain at the 21cm wavelength. The primary flux calibrators 3C 48 
(16.5 Jy), 3C 147 (22.5 Jy), and 3C 286 (15.0 Jy) were observed regularly to monitor the 
instrumental performance, and the antenna gain of 1.65±0.05 K/Jy is derived from these 
data. Our regular observation of bright calibration sources has verified the stability of the 
telescope gain factor. The rms noise after 60 minutes of on-source integration was 1.4 mK 
when the spectra were smoothed to 50 kHz (~10 km s _1 ) resolution. The corresponding 5a 
H I mass detection limit at a distance of 100 Mpc is 1.65xl0 8 100 M where AV is the 

line-width of the H I emission. 
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The GBT data also has very good stability within the bandwidth and the baselines 
could be fitted using first or second order polynomial for most of the HCGs. All the line-free 
channels with the exception of 200 channels at each end were used to fit baselines. Since 
most of the groups were observed near transit, we estimate an attenuation due to opacity of 
the air to be around 1% and corrected the GBT spectra accordingly. 

The second IF was tuned to search for the 1665 MHz and 1667 MHz OH maser transi- 
tions simultaneously with the H I transition. Dense molecular gas associated with nuclear 
starbursts or circum-nuclear disks h as be en detected in OH maser emission in nearby and 
distant galaxies (see review by iLol 120051 ) . Excited primarily through IR pumping, these 
transitions trace the presence of a substantial quantity of warm and dense molecular mate- 
rial (T = 40-200 K, n = 10 4-7 cm -3 ) heated by an intense radiation field. No OH maser was 
detected in any of the groups. Thi s result is consistent with the ab sence of infrared bright 
starburst systems in these groups (jVerdes— Montenegro et al.lll998l ). Thus, we concentrate 
on the analysis of the H I emission in the remainder of the paper. 



3. RESULTS 
3.1. GBT Data 

Among the 26 HCGs observed using the GBT, H I emission was detected in all but 
one group. The H I masses for the groups are calculated using the relationship M(HI) = 
2.36xl0 5 .D 2 (S'A\^) Mq where D is the luminosity distance to the group in megaparsec (Mpc) 
and the SAV is the velocity integrated H I flux density in Jy km s -1 . The derived H I masses 
within the GBT beam (~9.1') range between 0.5-28.5 x 10 9 M (see Table El for individual 
values including H I mass and linewidth). The velocity range over which the H I line flux is 
integrated is also presented in Table [21 The only undetected group, HCG 35, is the highest 
redshift group and is not a member of our distance-limited complete sample. The sensitivity 
of the GBT makes the 3cr upper limit for this high redshift group (cz = 16249 km s _1 ) 
M(HI) < 28 x 10 9 M Q assuming a line width of A V = 600 km s _1 (average line- width for 
our sample). We adopt a systematic uncertainty in calibration and baseline subtraction of 
10%, which is substantially larger than the measured noise level (around 1%) in most cases. 
For HCG 40, which was observed under the pilot program GBT05A-051, the antenna gain 
was not measured. We adopt the same gain factor from our observation under the program 
GBT05C-005. The reduced and calibrated spectra are shown in Figure [U 



Visual comparison of ou r GBT spec t ra wi th previous single dish measurements by 
Williams &: Rood! (Il987| ) and iHuchtmeierl (119971 ) confirm our calibration accuracy. A de- 
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tailed comparison is not discussed here. The large aperture, clean optics, and the low noise 
receiver of the GBT lead to superior calibration, sensitivity, and baseline stability compared 
with the earlier single dish measurements. This makes the GBT an excellent instrument for 
detecting the diffuse neutral component of the IGM, which may be very faint, extended, or 
broad in velocity (AV = 1000 km s _1 ). Also, the noise level ach ieved in our observ ations 
after an hour of integration is about 5 times lower than that of iHuchtmeierl (119971 ). The 
extended nature of the H I emission in these groups also makes a meaningful comparison 
with other single dish measurements with different beam sizes difficult. Instead, we focus 
our comparative analysis on the H I data obtained using the GBT and the VLA. 



3.2. Re-evaluation of Hi-deficiency for the Groups 

HI deficiency is a measure of the deviation of the observed H I mas s from the predicted 



H I m ass for a galaxy based on its optical luminosity and morphology (see lHaynes fc Giovanelli 



1984 Table 5, page 791). It is defined as Def HI = log[M(HI) pr . edicted ] - log[M(HI) o6ger „ ed ]. 
The predicted H I mass for a galaxy group is estimated as the sum of the predicted H I 
masses for individual galaxies (VM2001). 

The values of predicted H I masses used in ou r analysis are t aken f rom VM2001, which 



used optical luminosities and morphologies from iHickson et al.l (119891 ) . HCG 90 was not 
included in the VM2001 sample and thus we estimate the predicted H I mass for this group 
following the same procedure as VM2001. In most cases the GBT beam covered all the 
galaxies in a group, with the exception of HCG 10, 44, 68, & 93 where we correct the 
predicted H I mass to account for the GBT beam coverage. We assume that the predicted 
H I mass in individual galaxies is distributed uniformly inside the optical extent of the galaxy 
and hence the correction in the predicted mass translates to the fraction of the optical galaxy 
covered by the GBT beam. The observed H I masses, the predicted H I masses within the 
GBT beam and their corresponding H I deficiencies are provided in Table [2j 

A comparison of the predicted H I mass and the observed H I mass for our entire sample 
is shown in Figure [2j We find that most of the groups are H I deficient. While the predicted 
H I masses for most of the groups are similar, ranging from 9.34 to 10.57 dex; the observed 
H I masses vary over two orders of magnitudes due to the range in H I deficiency. The 
highest redshift group HCG 35 (part of our additional sources) can be seen with an arrow 
indicating the observed mass as a 3 sigma upper limit. H I was detected in the other three 
groups part of our additional sources and all of them lie in the H I normal region of the plot. 



The H I content and spatial distribution in the groups is dependent on the evolutionary 
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history of individual galaxies. Any transformation of neutral hydrogen as a result of tidal 
interaction within the group environment is expected to have some correlation with the 
group's H I deficiency. Therefore, we have divided the HCGs from our complete sample into 
three classes based on their predicted H I mass and the GBT observed H I mass: 



• normal H I content: observed mass more than 2/3 of the predicted value (HCG 23, 25, 
79, & 68); 

• slightly H I deficient: observed mass ranges between 2/3 to 1/3 of the predicted value 
(HCG 10, 15, 16, 31, 37, 40, 58, 88, 91, 92, 97, & 100); and 

• highly H I deficient: observed mass less than 1/3 of the predicted value (HCG 7, 30, 
44, 67, 90, & 93). 



These distinctions are based on the observed sc atter in the H I content for each Hubble 
type (~0.2 in dex; see lHaynes fc Giovaneililll984l ). This classification enables us to examine 
quantitatively the hypothesis that the H I content of each group indicates their age or 
evolutionary stage. 

The H I deficiencies estimated by VM2001 using previous single dish measurements are 
also presented in Table [2]for comparison. For 80% of the groups, the GBT detected more H I 
than previous single dish observation s using NRAO's 91m tel escope even though the beam 
size for these observations was 10.8' (IWilliams fc Rood! 119871 . VM2001). For three slightly 
deficient groups and two highly deficient groups the previous single dish values are larger 
than that measured using the GBT. The most likely source of this difference is the extended 
nature of the H I distribution which was captured better by the 10.8' beam as opposed to 
the 9.1' GBT beam. The histogram of H I deficiency evaluated from GBT data is shown 
in Figure [3j The peak of the distribution lies in the deficiency range of 0.2 to 0.4. The 
dotted curve represents the expected distribution of H I deficiency for a H I normal sample. 
This shows that the uncertainty in deficiency estimation cannot account for the observed 
distribution and therefore, confirms the existence of H I deficiency in HCGs. 

Even though these highly sensitive measurements decrease the H I deficiency of HCGs 
considerably, however, it does not elevate the total H I content of th ese groups to the n orma l 
level. The addition of hot IGM traced in X-ray may in some cases. iRasmussen et al.l (120081 ) 
analyzed new and archival X-ray data from Chandra and XMM- Newton telescopes for eight 
HCGs, and their X-ray results add a few times 10 10 to 10 11 M Q of hot IGM gas to the total 
gas content in some of the X-ray bright groups. This result suggests a direction for solving 
the H I deficiency problem by accounting for colder tidally stripped material that might have 
been ionized in the IGM. However, it is crucial to account for all sources of hot gas in galaxy 
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groups like accretion of gas from large-scale structures by the group potential, feedback from 
constitute galaxies as well as the deposition of cold ISM in the IGM as a result of tidal 
activities. It is also clear that the tidal debris is not the major contributor to the hot IGM 
gas mass as the missing H I mass (M( H t) wed - M( H l) & s ) for X-ray bright groups like 
HCG 15, 37, 40, fc 97 is on an average an order of magnitude smaller that the hot gas mass 
( iRasmussen et al.ll2008l ). 



3.3. Comparison with the VLA Data 



It is a common practice to compare the flux recovered by an interferometer with a single 
dish measurement in order to characterize the limitations of the resulting aperture synthe- 
sis images. For instance, an interferometer such as the VLA is a spatial filter and is thus 
insensitive to diffuse, extended structures > 15' in the D-arrajfl . Although the collecting 
area of the GBT and the VLA are comparable, the VLA has a much lower surface bright- 
ness sensitivity than the GBT. In addition, the GBT has a larger instantaneous bandwidth 
(A V ~ 2500 km s _1 ) that allows identification of emission features with a much broader 
(> 500 km s _1 ) line-width. In comparison, the bandwidth of the VLA is less than half 
of the GBT, and broad spectral features can be easily missed by the VLA due to limited 
spectral coverage. 

The VLA data used for comparison come primarily from our own imaging survey 



of HCGs ( 


Verdes— Monteneero et al. 


2001; 


Williams et al. 


2002 


2005; 


Borthakur et al. in prep 


2010; 


Yun et al. in prep 


2010). 



for six groups are constructed by our own reduction of the archival data that were pub 
lished previously (IWilliams fc Van Gorkoml 11995c IWilliamsl Il998l ; iBarnes fc Websterl 12001 



Heiderman et al.l 120051 ) . Table |3] summarizes the details of the VLA data used for compar- 
ison. The flux density sensitivity of the comparison VLA data is similar in most cases (~ 
0.5 mJy/beam), but the angular and spectral resolutions vary slightly. 

Comparisons of the new GBT spectra with the spectra derived from the VLA obser- 
vations for our complete sample as well as for two of the additional sources are shown in 
Figure HJ We define the difference between the VLA and the GBT spectra as the excess H I 



For a uniformly weighted, untapered map produced from a full 12 hour synthesis observation of a source 
which passes near the zenith, the largest angular scale that could be detected by the VLA in D-array is 15'. 
However, since most of the observations are shorter than 12 hours and used natural weight the maximum 
angular scale should be smaller than 15' 
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mass 



M, 



excess 



M G bt — M VLA 



(1) 



We make a quantitative characterization of the excess gas as a fraction of the total H I in a 
group and define it as the excess gas mass fraction, / 



In groups where the excess H I mass dominates, the value of this ratio is greater than 50%. 
Similarly, in groups where this ratio is smaller than 50%, the primarily H I component is the 
high surface brightness H I detected by the VLA. 

Figure [5] is a histogram of the total H I mass and the excess H I mass in our sample. 
The excess gas mass fraction,/ for our complete sample varies between 5% and 81% with an 
average value of 36% (Table HJ). In order to ensure that this comparison is fair, we multiply 
the original VLA map with the GBT beam pattern and then summing the H I flux in each 
channel to construct the corresponding VLA spectra. Therefore, the VLA spectra used in 
the comparison are for regions identical to those observed by the GBT. Furthermore, the 
quantitative comparison of the H I mass is limited to the velocity range that is common 
to both spectra although the GBT detected much broader emission features in some cases. 
Consequently, the values mentioned above are lower limits of the H I missed by the VLA. 
The GBT data is also much more sensitive than the VLA data and this easily noticeable in 
HCGs with faint H I emission such as HCG 15, 37, 90, and 97. 

Extra steps were taken to ensure consistency in the calibration as reliable cross-calibration 
is critical for quantitative comparisons. To ensure excellent calibration with respect to the 
VLA, we observed the same three VLA primary flux calibrators 3C 48, 3C 147, and 3C 286 
throughout the observing run. By taking these extra calibration steps, we estimate that the 
calibration between the GBT and the VLA spectra should agree within 10%. The excellent 
agreement between the GBT and the VLA spectra for groups dominated by high surface 
brightness (HSB) features, such as HCG 7, 10, 16, and 67, nicely demonstrates the success 
of our cross- calibration. Even in groups where there is an overall difference in total flux such 
as HCG 26 and 79, parts of the spectra match well within the statistical errors. The fact that 
the total flux measured by the GBT always agrees or exceeds the VLA measurement (except 
in the case of HCG 40 where we have a large calibration uncertainty) but is never smaller 
rules out the possibility of relative calibration accuracy varying wildly. Therefore, calibration 
errors are not responsible for the dramatic differences seen in HCG 31, 44, and 100, and the 
excess emission in the GBT spectra is real. The excess emission is also correlated with the 
group properties (see § I4.4p . and this is a strong indication that the measured excesses are 
real. One may plausibly suspect a continuum subtraction error for missing broad spectral 



/ 



M GB T — My LA 
McBT 



M, 



Mqbt 



excess 



(2) 
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features in the VLA data. However, the absence of negative flux values in the VLA images, 
which is expected from overdoing continuum subtraction, rules out significant errors in VLA 
data reduction. 

Three broad categories on the spectral distribution of the excess gas have emerged from 
our comparison of the GBT and the VLA spectra: (I) identical line widths but with a large 
line flux difference; (II) broad emission wings (50-200 km s -1 ) associated with VLA detected 
spectral features; and (III) broad (> 400 km s -1 ) emission features outside the VLA H I 
line-width. The most commonly seen is the category I, and nearly all groups show this trait 
to a degree. The identical line widths is a strong indication that the additional emission 
detected by the GBT is associated with the structure traced by the VLA. It is possible that 
additional, independent H I systems sharing the same velocity range are present, and a few 
distinct features with a narrow line-width were seen in the GBT spectra. 

The second category with broad emission wings is also relatively common. HCG 92 
("Stefan's Quintet") and HCG 100 represent two outstanding examples of this type, charac- 
terized by broad H I emission wings in excess over the velocity range detected by the VLA. 
The wing features in HCG 92 cover the velocity range 6300 to 6530 km s _1 and 5750 to 
5900 km s _1 where as in HCG 100 the wing feature can be seen in the velocity range from 
5450 to 5650 km s _1 . The excess can be seen in the GBT data even when it is binned to 
match the resolution of the VLA spectra (see Figure [6]). A collision involvin g an intruder at 



a vel ocity in excess of 1000 km s 1 is currently occurring in HCG 92 (see lAppleton et al. 



20061 . and references therein), and the broad H I wings are plausibly associated with high 
velocity shocks in this case. The VLA imaging of HCG 100 shows a cloud of H I surrounding 
most of the member galaxies and extending far beyond the group diameter in the VLA image 



(IBorthakur et al. in prep 1 120101 ) . The VLA spectrum is clearly missing the broad wings on 



both sides of the main H I peak. Other groups showing similar broad emission wings are 



HCG 31, 44, and 79. 



The third category is the broad, faint emission features unseen in the VLA spectra. 
In HCG 15, we detected 56% of total H I in this system as excess H I emission spread 
over a velocity range of more than 1000 km s~ l . The VLA channel maps of H CG 15 show 



H I emission associated with a single galaxy (IBorthakur et al. in prep 1120101 ) spread over 
no more than 250 km s _1 . The VLA not only missed the broad H I emission because of 
a smaller bandwidth, but also could not detect fainter emission due to the limitation in its 
surface brightness sensitivity. The VLA o bservations of HCG 30 an d 37 have failed to detect 



any high surface brightness H I emission ( lYun et al. in prep 1120101 ) while their GBT spectra 



show faint, broad emission covering the entire redshift range of the member galaxies (see 
Figure H]). These observations suggest that the entire H I emission in these groups is in a 
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diffuse component. In addition, HCG 44, 58, 79, 90, and 97 also show evidence of faint, 
broad emission not detected by the VLA. 

HCG 31 is one of the most interesting cases where the comparison of H I spectra 
shows an enormous difference between the GBT and the VLA, confirming the presence of 
a substantial diffuse neutral IGM with much larger H I mass than was previously believed. 
The peak flu x and the line pr o file o f the GBT spectrum matches very well with spectra 
published by [Williams fc Rood! (119871 ). thus confirming our calibration. Another interesting 
fact that stands out in the comparison plot for HCG 31 in Figure H] is that discrete features 
seen in the VLA spectrum are not seen in the GBT spectrum. In fact, the GBT spectrum 
is much smoother, and the velocity range between the discrete peaks at 4000 km s _1 and 
4120 km s^ 1 is filled in. It is likely that the diffuse H I emission spread over a velocity 
range AV ~ 200 km s -1 accounts for the majority of the GBT spectrum. The diffuse H I 
is likely associated with features much larger in angular extent than group size, possibly in 
H I filaments being accreted by the group, simi lar to those seen in the VLA image reported 
previously by I Verdes— Montenegro et al.l (120051 ). Additional observations are underway for 
an in-depth investigation of the nature of the excess H I detected by GBT. 



4. DISCUSSION 



H Icontent in HCGs varies enormously from one group to another. VM2001 found a 
broad range of H I deficiency with a mean H I content of 24% of the expected amount 
(Def H i = 0.62 ± 0.09) for the individual galaxies and about 40% (Def m = 0.40 ± 0.07) for 
the groups. The H I distribution in these groups is also quite diverse, ranging from groups 
dominated by normal galaxies with distinct individual H I disks to groups with galaxies 
submerged in a single common H I cloud. Evidence for strong tidal interactions is frequently 
seen in the H I distribution imaged by the VLA. Such interactions continuously shape the 
gas distribution and morphology . Consequently, the H I content and mo rphology reflects the 
evolutionary history of a group (1 Verdes— Montenegro et al.ll2000l . 120011 ). 



A new result from our new GBT survey is that the excess H I emission is quite common 
among HCGs. Here we attempt to understand the distribution and nature of the excess gas 
through a quantitative analysis of the VLA and the GBT spectra and group properties. 
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4.1. Geometry and Kinematics of the Intragroup HI 

The amount of H I detected by the GBT exceeds the VLA detected amount by a factor 
of two or more in HCG 15, 31, 37, 44, 68, 92, 97, and 100, suggesting that the majority 
of H I in these groups is associated with extended diffuse structure, undetected by the 
VLA. We experimented with the VLA data to look for signs of a spatially extended H I 
component using multi-scale CLEAN in the IMAGR procedure of the Astronomical Image 
Processing Softwa re (AIPS) which enhances flux associated with larger angular scales (or 



shorter baselines) (jCornwellll2008l ; iRich et al.ll2008l ). Figure [7] shows an overlay of the multi 



scale CLEANed VLA spectra of HCG 100 over its GBT and the VLA spectra. The VLA 
multi-scale CLEAN spectrum was obtained by using three Gaussian source model widths- 
0, 200 and 300 arcsecs, in deconvolving the images. Even in this case, not all of the GBT 
flux is recovered, and nearly all "recovered" fe atures are faint extensions of fea tures already 



seen in the standard reduction map (details in lBorthakur et al. in prep II2010I ) 



Thus, the missing gas in the VLA images may exist as either: (1) low column density 
diffuse medium which can only be detected by GBT due to its superb surface brightness 
sensitivity; (2) spatially extended sheets, filaments or extended disks that are missed by the 
VLA because of spatial filtering; or (3) as a combination of both. A physical parameter that 
provides a quantitative way of distinguishing between these paradigms is the filling factor 
of the excess gas. We characterize the filling factor as the fraction of the GBT beam the 
excess gas would cover if it was distributed uniformly with column density corresponding to 
the noise level (3 ctn(hi),vla) in VLA maps. Thus, the filling factor A// is computed as 

-^-excess -^excess 

^ 3 <JN(HI),VLA &GBT Mh M SaV LA 

where floBT is the area of the GBT beam in cm 2 and M# is the mass of an hydrogen atom. 

Table H] summarizes the excess H I masses, VLA noise properties and the filling factors 
for each of the groups. The 3a column density (seventh column of table H]) was derived 
assuming that the features are spread over 2 channels thus implying an intrinsic line-width 
of 42 km s -1 . The derived filling factor covers a wide range of values, from as little as 3% 
in HCG 79 to 566% in HCG 68. For most of the HCGs, A ff is much smaller than 100%, 
which supports the presence of a faint diffuse H I medium that could have been missed by 
the VLA due to its lower surface brightness sensitivity. For groups with Aff > 100%, spatial 
filtering by the VLA leading to a significant loss of flux is a highly plausible explanation. 
This is corroborated by the fact that in all four cases where the filling factor is greater than 
100%, (HCG 44, 68, 91, & 100), VLA maps show cloud-like structures or extended disks 
covering around 8-10 ' (nearly the entire GBT beam). Furthermore, striking similarities in 
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the spectral profile between the VLA and the GBT spectra in these five systems strongly 
advocate that the excess gas is dynamically similar to the H I structures mapped by the 
VLA. 

The excess gas is often associated with a broad; smooth spectral feature rather than 
discrete peaks. The dominant spectral feature in spectra of H I deficient groups like HCG 15, 
90, and 97 is a faint broad feature. The GBT also detected broad wing-like features in many 
other groups including HCG 92, 79, and 58. A noteworthy characteristic is that these broad 
emission tends to fill the gaps between peaks in the spectrum, resulting in a much smoother 
spectrum such as in HCG 31 (also see HCG 26, 67, & 100). This suggests the presence of a 
diffuse fainter H I component distributed over a large velocity range. The probable origin of 
this H I component is the evolution of tidally stripped ISM in the hot IGM (see § I4.3p . A 
detailed discussion of the fate of tidally stripped ISM clouds in the hot IGM is presented in 
the next section. 

In summary, we conclude that the excess gas is spread over a large velocity range and 
is smoothly distributed as a diffuse, low surface brightness medium. The broad spectral 
signatures specifically rule out cold discrete filaments. The filling factor values suggest that 
the excess gas is extended and must occupy at-least a significant fraction (median = 40% ) 
of the GBT beam. Consequently, the excess gas exists as a faint diffuse neutral intragroup 
medium in most HCGs. In a few groups where Aff > 1, large HSB H I structures that 
were spatially resolved out by the VLA may also contributes to the excess. A combination 
of the diffuse IGM and some HSB features is likely in groups exhibiting large line-flux 
differences along with broad wing-like spectral features (e.g., HCG 100). These groups may 
also represent a transitional phase occurring between groups with clumps of tidally stripped 
H I in their IGM and groups with a faint diffuse extended H I component distributed in their 
IGM. 



4.2. Lifetime for H I Structures in IGM 



Compact groups are evolving structures. The dynamical time {td yn ) for a typical group 
with a diameter of 100 kpc and velocity dispersion of 250 km s _1 is around 400 Myr. During 
this time, galaxies experience multiple tidal interactions with each other, and a considerable 
amount of ISM may be removed from individual galaxies and deposited in the group envi- 
ronment. This is evident in groups like HCG 16 which show multiple tidal tails (VM2001) 
and also in groups like HCG 92 where the H I has been tidally removed from their parent 
galaxies and deposited as large H I clouds in the IGM (jWilliams et al.ll2002l ). The nature 
and fate of neutral gas in such extreme environments are determined by the interplay of 
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heating and cooling in these systems. The fate of tidally stripped debris can provide clues 
about possible causes of H I deficiency commonly observed among galaxies in such dense 
environments. 

In this section, we attempt to check theoretically if the neutral tidal debris should survive 
for a considerable amount of time for us to be able to detect it in systems like the HCGs. In 
other words, should we expect to find a surviving component of the tidal debris deposited in 
the IGM? In order to do so, we investigate the timescales and cond itions for the survival o f 



neutral clouds in the hot IGM. For a spiral dominated galaxy group, iMulchaey et al.l (jl996bl ) 
calculated the expected virial temperature of the IGM to be around 2 x 10 6 K. Any cold H I 
cloud/filament will be conductively heated by the hot IGM. This can ionize the cloud and 
incorporate it into the hot IGM. In addition, a neutral cloud will expand into the IGM due 
to its intrinsic velocity dispersion. This will eventually lower its column density and make 
it more susceptible to ionization by the ultra-violet radiation background. Based on the 
timescales of these processes, we can derive the expected lifetime of a typical tidally stripped 
cloud. Since we have ample evidence of diffuse neutral IGM in HCGs, we can explore the 
parameter space for physical conditions that will ensure the survival of the H I clouds over 
the dynamical time-scale of the groups. This will help in constraining the properties of the 
H I clouds as well as the IGM, which in most galaxy groups escape detection by being too 
cold to emit in the X-ray. 



4-2.1. Conductive Heating Lifetime 



We begin by investigating the effects of conductive heating on a cold neutral cloud 
embedded in hot IGM. This process is not instantaneous, and the evaporation timescale 
depends on the physical properties of both the H I cloud as well as the IGM. Conductive 
heating can take place in two paradigms, classical or saturated, depending on the ratio 
of the mean free path of the electrons to the temperature scale-height of the hot IGM 
(jCowie &: McKed 119771 ). For a spherical neutral cloud with radius R c ; (in parsecs) and 
density n c ; (in cm -3 ) floating in the hot IGM at a temperature of Tjgm (in Kelvin) and 
density utcm (in cm " 3 ), the ratio of classical heat flux to saturated heat flux is given as 



flVollmer et al.ll200ll eq. 11 ) and , 



o-o ~ (T /GM /1.5 x 10 7 ) 2 /(n IGM x R d ) (4) 

The classical evaporation is applicable for systems where a a <C 1, where as saturated 
evaporation will be applicable for systems with a a ^> 1. The evaporation timescale in the 
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classical paradigm is given by (jCowie &: McKeelll977l eq. 22) 



t c ev = 3.3 x 10 20 n d Rl T^g 



Ink 
l30~ 



yrs 



(5) 



where In A = 29.7 + In n J( i,^(T/ GA //10 6 K). For a fixed IGM temperature and cloud mass, 



Kv °c n d Rl oc 



Mc, 

R r 



(6) 



This implies that a smaller and consequently denser cloud will much survive longer than a 
less dense cloud of the same mass. 



In the saturated paradigm, the evaporation timescale is given by (jCowie McKeeHl977l . 
eq. 64) 

Kv ~ !0 6 n cl nJ GM R d yrs (7) 

From Eq. [5] and [7J it is clear that the lifetime of the cloud increases as its density and/or 
radius increases in both paradigms, although at different rates. Also the lifetime is inversely 
related to the temperature of the IGM. 

Many of these physical quantities are poorly understood for compact groups, but we can 
make educated guesses based on our understanding of tidal interactions, I SM phases, and 



the in tergalactic medium. From multi-galaxy interactions seen in HCG 92, 1 Williams et al. 



(120021 ) found H I structure associated with HCG 92A to be as small as 7.3 x 4.4 kpc 2 . Thus 
for our calculations we assume a lower limit of tidally stripped neutr al cloud to be around 
1 kpc. In a recent study of jet-IGM interactions, Freeland et al] (2008) derived IGM density 
between 9 x 10~ 4 and 4 x 10~ 3 cm -3 . Using X-ray data, iRasmussen et al.l ( 120081 ) derived 
IGM densities for HCGs that ranged from 6 x 10~ 5 to nearly 4 x 10 -2 cm -3 . Adopting 
n IGM =10- 3 cm" 3 , T IGM = 2 x 10 6 K, and R d 
classical paradigm is applicable. 



10~ 3 cm" 3 , T T hm= 2 x 10 6 K, and R d = 1 kpc, we find a = 0.0177 < 1 and the 



Tidally stripped H I clouds have been part of the ISM of the individual galaxies and 
had similar de nsities and ionization fract ions as the H I emitting clouds typically found in 
the ISM. From iMcKee fc Ostrikerl (Il977l ). the lowest density phase of ISM with ionization 
fraction significantly less than 1 is the warm ionized medium (WIM). The density and the 
ionization fraction in this phase is 0.25 cm -3 and 0.68 respectively. For phases with lower 
ionization fraction than 0.68, the density of the ISM will be much higher. Using H I cloud 
density of 0.25 cm -3 and the above-mentioned values for the ISM and IGM parameters, we 
derive an evaporation timescale of ~ 15 Gyrs, which is 3> td yn . In the classical paradigm, 
the evaporation timescale is strongly correlated with the size of the cloud. The lower limit 
of H I cloud size that can survive the IGM for more than td yn is ~200 pc. For instance, the 
evaporation timescale for a smaller cloud with radius 100 pc is ~ 0.15 Gyrs. 
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Constraining the physical parameters of the IGM 

Our GBT detections, combined with information from VLA images, confirm that H I in 
HCGs survives long enough to be seen frequently in their harsh environments. Using this 
crucial piece of information, we can constraint the physical conditions of the H I clouds as 
well as the hot IGM by finding the region in the parameter space where the evaporation 
timescale is greater than the dynamical timescale of the groups. Throughout our analysis 
we will assumed the H I clouds to be stationary. If the proper motion of the H I clouds could 
be estimated more stronger constraints could be achieved. Unfortunately with the present 
state of observational data such estimation is not possible. 

Figure [8] shows the cloud evaporation time scale as a function of density of the IGM 
and the radius of the cloud. The solid line divides the parameter space into the classical 
and saturated evaporation regimes. The white region represents the range for which the 
evaporation timescale is more than 400 Myrs (typical dynamical time) where as the grey 
region is the range of parameters for which a neutral cloud would not survive the heating 
and evaporate within a single crossing time. The timescale increases or decreases as the 
parameter values moves away from the solid line and follows the proportionality seen in 
equation [H] and [7J It is evident that for clouds of radius >200 pc with an average density > 
0.25 cmT 3 are robust against evaporation for an IGM with T IGM = 2x10 s K, irrespective of 
niGM ■ 

The large timescale for heating can also be understood in the context of the critical 
radius, which is defined as the cloud size at which radiative cooling balances heating by 
thermal conduct i on. F or an IGM with Tjgm ^ 10 5 K, the critical radius as formulated by 



McKee fc Cowid (119771 . eq. 14 ) is 

R c « 4.8 x 10 5 Tf GM /n IGM cm. (8) 

The critical radius for a cold cloud in a hot (Tigm — 2 x 10 6 K) medium is 0.6 ^ R c ^ 6 kpc 
when the density of the IGM is 10 -3 Uigm ^ 10 -4 cm -3 . Thus, large H I structures with 
Rd > R c should survive for more than td yn (400 Myrs) due to a balance between conductive 
heating and radiative cooling. 

We have also explored the dependence of t ev for a small range of IGM temperatures and 
HI cloud densities. Figure [9] shows the parameter for eight different combinations of IGM 
temperature and cloud density. The plots have been arranged with increasing temperature 
from top to the bottom and increasing density from left to right. These plots clearly show 
that evaporation timescales decrease with increase in the IGM temperature and with the 
decrease in the density of the cloud, as expected from Eqs [5] and [71 Given that we detected 
21 cm emitting cold gas in most of the groups, we can conclude that the IGM temperature 
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is less than 4 x 10 6 K unless the density of the hot IGM is unusually low. Also, the density 
of the H I clouds, responsible for excess emission in the GBT spectra, is most likely higher 
than 0.01 cm -3 unless the IGM is substantially cooler than 2 x 10 6 K. This follows from 
the fact that the clouds with density lower than 0.01 cm -3 would not survive long enough to 
be detected in an IGM corresponding to typical virial temperature of 2 x 10 6 K. Since the 
relevant parameter space is mostly in the classical regime, the variation of IGM temperature 
has a much stronger effect on the timescale than the density of the cloud. On the other 
hand, the density of the IGM effects the timescale only in the saturated evaporation regime 
and consequently, we cannot draw a strong constraint on the yiigm from our data. 



4-2.2. Lifetime Based on Expansion of the HI structures 



In addition to heating, the tidal H I debris should also experience expansion as a result 
of their intrinsic velocity dispersion unless they are bound by their self-gravity. Such an 
expansion of an H I cloud would lower its column density, making it prone to ionizatio n 
by the intergalactic UV r adiation field. The 21 cm H I studies by ICorbelli et all ( 119891 ); 
Corbelli & Salpeterl (119931 ). and others have shown that a minimum column density of 2 x 



10 19 cm -2 is required for an H I cloud to shield itself from the background UV radiation 
field. Therefore, the expansion timescale, which refers to the time an H I cloud would take 
to reach the above mentioned column density, is important in estimating the lifetime of an 
expanding tidal structure. As a cloud expands its column density decreases which can be 
expressed in terms of the initial column density (JVj) and initial radius Ri as 



N 



Area; 



{Ri 



Area (Ri 



+ v t e ^ 2 



(9) 



Assuming that the tidally stripped cloud in question is a disk/cylinder with Ri =1 kpc 
and Ni = 3 x 10 21 cm -2 , expanding at v = 20 km s . The expansion timescale for this 
cloud to reach the above mentioned threshold value is around 500 Myrs. In this time the 
cloud expands more than 100 times its initial area (or 1000 times its initial volume), thereby 
creating a faint diffuse neutral IGM. 



4.3. Possible Origins of the Diffuse H I Emission 



Based on the observed deficiencies as well as complex H I tidal structures seen in VLA 
imaging of the HCGs, the excess gas mass can be attributed to evolving tidal debris deposited 



- 19 - 



by the individual galaxies. The fact that the excess gas does not render the HCGs super- 
abundant in H I supports the tidal debris origin of the excess gas. Furthermore, the velocity 
distribution of the excess gas is extremely similar in line shape to that of the VLA detected 
high surface brightness H I associated with individual galaxies (catagory I) in more than 80% 
of the groups. The excess gas content also increases with the evolutionary phase (VM2001) 
of the group (discussed in detail in § I4.4.3[) . which supports the idea of evolution of tida l 
deb ris into this component. Similar detections of tidal H I structures by iDavies et al.l (120041 ) 
and lOsterloo fc van Gorkoml (120051 ) in an even more extreme environment of the Virgo inter- 
cluster medium also lends support the existence and preservation of tidally stripped neutral 
gas for > 10 8 years. 

A potentially important process for the IGM in some of the X-ray bright compact groups 
such as HCG 97 and HCG 15 is the preci pitation of col d gas from the hot IGM through the 
radiative cooling (so-called "cooling flow" IFabianlll994l ). Unless a constant source of energy 
is present to reheat the gas back to the IGM temperature, hot gas is expected to cool down 
i n the centers of groups and clusters. Calculations based on work by iNarayan fc Medvedev 
( 1200 ll ) suggest that the transfer of heat from hot IGM to the condensing cores would take 
about 0.5 Gyr for a group of size 100 kpc. Such a conduction time-scale for the IGM 
suggests the possibility of a neutral IGM resulting from cooling of the X-ray bright IGM at 
the centers of these groups. However, most of the HCGs used in our study are not X-ray 
bright and therefore, cooling flow is not likely to be the dominant mechanism responsible for 
the existence of the neutral gas in our sample. Therefore, we can safely conclude that the 
majority of the excess gas detected by the GBT is not a result of condensation of the hot 
IGM. 



In a handful of groups where the GBT beam covered a region much large r than the group 
size, accretion of cold gas from the large-scale structure (IKeres et al.l 120051 ) is a possibility. 
For most of the HCGs we have surveyed the central region of the groups where the accreted 
cold gas, if any, is expect to get shock heated and ionized thus eliminating the possibility of 
its survival for ubiquitous detection in 21 cm emission. Therefore, the most plausible origin 
of the excess gas found in the central regions of the HCGs is H I tidal debris deposited in 
the IGM as individual galaxies interacted with each other and the group potential. In this 
scenario, no additional cold gas contribution from accretion or cooling flows (for X-ray bright 
HCGs) is required, although minor contributions from these processes cannot be ruled out. 
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4.4. Excess H I and Group Properties 

^-excess is substantial in most of our groups and has important implications on our 
understanding of the group medium and its evolution. As discussed earlier, this excess gas 
most likely exists as an extended diffuse medium that is being created and maintained by 
depositing tidally stripped ISM from individual galaxies into the IGM. In this section, we 
explore the distribution of the excess component using excess gas mass fraction in terms of 
the group radius, the H I deficiency, the evolutionary stage as well as Hubble type of the 
central galaxy of the HCGs. For our statistical analysis, we use a sub-sample of our complete 
sample for which VLA D-array H I data was available for comparison. The additional sources 
- HCG 26 and HCG 48 are not included. HCG 40, although a part of our complete sample, 
is also excluded from the analysis due to the uncertainty in the calibration of the GBT data. 
In addition, VLA spectra has more flux than in the GBT spectra for HCG 40 and following 
equation [21 fnccm < 0. Hence, our statistical analysis sample consists of - HCG 7, 10, 15, 
16, 23, 25, 30, 31, 37, 44, 58, 67, 68, 79, 88, 90, 91, 92, 93, 97, and 100. 

4-4-1 . Radial Dependence with Group Diameter 

HCGs represent a diverse set of galaxy groups with a range of sizes, number of members, 
H I content and evolutionary phases. In our sample, the angular diameter of the groups varied 
between 1.3' to 16.4', which translates to 0.14 to 1.80 times the GBT beam at 1.4 GHz. The 
excess gas mass fraction is shown as a function of the group diameter sampled by the GBT 
beam in Figure (TUl For groups where / > 0.5, the GBT beam covers the inner region up 
to a maximum of 3 times the group diameter. For instance in HCG 44, where the excess 
gas mass fraction is 60%, the GBT beam covers only the inner one-third of the group area 
(group diameter is 16.4'). In other words, the diffuse excess gas is the major contributor to 
the total H I mass in the central region of this system. This indicates that the distribution 
of the diffuse component is within the central region of the groups along with the individual 
galaxies. The majority of groups are either lack any diffuse component or are dominated 
by it with few groups lying in between as seen in Figure [10J In fact, this offers the first 
indication that the / distribution is clustered, with one peak near / ~ 0.2 and a second 
peak near / ~ 0.6, irrespective of the fraction of group diameter covered by the GBT. There 
are only two systems with 0.3 < / < 0.5. This suggests that the transition from the low 
diffuse gas phase to high diffuse gas is not gradual. Instead, the transition may be quite 
rapid and the change in the ratio indicates an evolution in the group property. 

In terms of the physical dimension, the GBT beam coverage varies from 40 kpc to 
250 kpc. Figure [TT1 shows the same data as Figure [TOl but in physical units (kilo parsec). 
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The trends seen are essentially identical to those presented in Figure [TOj but a potentially 
important and new insight is offered: groups where the surveyed region is smaller than 100 kpc 
show a rather significant fraction of excess gas, f > 0.4. While the statistics are poor, such 
a trend is expected only in scenarios where the group radius actually reflects the physical 
proximity of the member galaxies. Specifically, such a trend should not exist in apparent 
groups identified by chance projection only (e.g iHernquist et al.lll995h . 



4-4-2- Dependence on El-deficiency 

Our HCG sample can be divided into three categories on the basis of their H I defi- 
ciencies, as normal H I content, slightly H I deficient, and highly H I deficient (see § 13. 2j) . 
The original definition of H I deficiency was based on the H I content of individual galaxies. 
However, in HCGs the H I distribution varies enormously from the distribution of the opti- 
cal galaxies, and it is often not possible to associate H I structures with individual galaxies. 
Therefore, we extend the definition of H I deficiency to groups by evaluating the ratio of the 
sum of the predicted H I mass for individual members to the observed H I mass associated 
with the entire group. 

The excess H I fraction for the three deficiency classes of HCGs is summarized in 
Table O The mean (median) excess gas mass fraction, /, for the normal H I content, slightly 
H I deficient, and highly H I deficient groups are 30 % (26%), 42 % (51%), and 28 % (20%) 
respectively. The comparison of H I deficiency for the groups is also shown graphically on 
the top panel of Figure [12j Most of the groups with a high fraction of excess gas are slightly 
H I deficient groups whereas most of the normal and highly deficient groups have much 
smaller fraction of excess gas. The scatter within each class is large which makes the average 
distribution ineffective in deducing any statistically significant conclusion. 

On the other hand, for groups with H I deficiency < 0.4 (= 2 o in predicted H I mass), 
one would expect the two tidally stripped H I components - diffuse H I detected only by 
the GBT and the high surface brightness (HSB) H I detected by the VLA, to constitute the 
total H I in these systems. Figure [13] shows a comparison of the excess gas content and the 
HSB H I content in terms of the predicted mass for a sub sample of 11 HCGs, where H I 
deficiency < 0.4. In general, the fraction of diffuse excess gas increases with the increase in 
the H I deficiency. The Spearman's (rho) rank correlation coefficient for this subsample is 
—0.523919 with significance of its deviation from zero of 0.098. This provides an indication 
that HSB structures in the IGM are evolving into the diffuse H I component before being 
ionized and assimilated into the hot IGM. 
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In terms of the evolution and fate of the neutral tidal debris, the diffuse gas can be 
understood as an intermediate stage between HSB H I and ionized hydrogen. In this picture, 
the stripped H I from the galaxies begin their life as tidal features (VM2001) in the IGM, 
which evolves into a diffuse faint cold gas component that was detected with the GBT. After 
a certain threshold, depending on the physical parameters of the IGM and the diffuse H I 
clouds (see § 14.2)) . the hot IGM ionizes and assimilates the diffuse component. Thus, in 
the early stages of evolution, groups show an increasing fraction of diffuse gas component 
with increasing deviation from H I normalcy. Towards the later stages of evolution, groups 
continue to show strong H I deficiency as the diffuse H I component is transformed into 
the ionized IGM and cold neutral gas is truly lost from these systems. Since the tidal 
interactions do not continuously deposit constant amount of neutral debris in the IGM, 
no correlation between the total HSB component and the total diffuse gas component is 
expected or observed for highly evolved systems. 



4-4- 3. Dependence on Group Evolutionary Phase 

HCGs represent the full range of formative and evolutionary stages of galaxy groups. 
VM2001 proposed an evolutionary sequence in HCGs based on their H I morphology and 
distribution. According to this model, during the first phase of evolution (Phase 1), the H I 
distribution and kinematics are relatively unperturbed, and more than 90% of the H I mass 
is found in the disks of individual galaxies with the remaining gas found in incipient tidal 
tails. In the second stage of evolution (Phase 2), groups still retain a significant amount of 
H I in the disks, but 30% to 60% of the total H I mass forms tidal features. In the final 
phase, most of the H I is stripped from the disks and is found in tails produced by tidal 
interactions or not detected at all (Phase 3a). A slightly less common phase is the Phase 3b 
where the H I forms a single large cloud containing all galaxies with a continuous velocity 
gradient and a single peak line profile. 

A statistical comparison of the excess gas mass fraction, /, for our sample based on 
their evolutionary stage is shown in Table [61 For this analysis we do not distinguish between 
Phases 3a and 3b. The mean (median) excess gas mass fraction for groups in the evolu- 
tionary phases 1, 2, and 3 are 21 % (21 %), 30 % (17 %), and 48 % (54 %), respectively. 
Dispersion in excess H I fractionamong the members in a phase may be understood in terms 
of a continuous evolutionary sequence in HCGs where the three phases are broadly (and 
somewhat arbitrarily) divided to quantify the gradual change. It is apparent that in the 
initial phases of evolution there is much less diffuse neutral IGM than in the last phase, 
where the neutral IGM is well developed. For instance, all the groups in the Phase 3, with 
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the exception of HCG 93, have a considerably large fraction of diffuse H I . HCG 93 may 
be a special case where the addition of a new H I rich member (93B) reduces the fractional 
significance of the diffuse neutral H I in the IGM. 

A graphical comparison of excess gas mass fraction shown in the middle panel of Fig- 
ure [12] demonstrates the clustering of Phases 1 and 2 towards lower values and Phase 3 
HCGs towards higher values of the excess gas mass fraction. As discussed in § 14.4. H this 
clustering may indicate a rapid transition for groups dominated by H I mostly in the ISM 
of individual galaxies to a phase where most of the H I is found in the IGM. We emphasize 
that the group evolutionary phase is determined from the H I morphology as mapped by 
the VLA while the excess gas mass fraction comes from the integral total line flux in the 
GBT spectra. One should not a priori expect any correlation between these two completely 
independent parameters, except for an underlying physical mechanism that governs them 
both. Groups with / < 40% in Figure [12] consist mainly of HCGs in Phases 1 & 2 while 
those with/ > 40%) are almost entirely from Phase 3. In the group evolutionary scenario, 
this trend indicates that the evolutionary Phases 1 & 2 are much closer in timescale and 
in terms of the neutral IGM build-up process, and the evolutionary phase 3 is indeed more 
distinct. Again, the small sample size and a large dispersion within each sub-sample is an 
important limitation. 



4-4-4- Central Dominant Galaxy and Diffuse Intragroup HI 



In the cold dark matter (CDM) paradigm, dark matter halo properties such as mass 
and angular momentum determine the galaxy properties. Specifically, for a given massive 
halo, the nature of the central galaxy is thought to be determined by the merger history 
of the halo. By examining the de pendence of the c ross-correlation between galaxies and 
galaxy groups on group properties, lYang et al.l (120061 ) found that the dark matter halo bias 



of galaxy groups decreases as the star formation rate of the central galaxy increases. These 
authors interpret this finding in terms of the halo bias dependence on halo formation time 
( iGao et al.l 120051 ) , in the sense that halos that were assembled earlier are more strongly 
biased with a redder central galaxy. In other words, on an average galaxy groups with a red 
central dominant galaxy have formed much earlier than groups with a blue central dominant 
galaxy even if their present halo masses are the same. 



An e xcess of early type galaxies in HCGs have been reported by iHickson et al.l ( 119881 ). 



Sulentid ( 120001 ) and VM2001 suggested that some of the lenticulars and SO galaxies were 
spirals that were stripped of their H I disks as a result of galaxy interactions. We have 
examined whether the excess H I content has any dependence on the Hubble type of the 
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group central dominant galaxy. If the earlier Hubble type of the central dominant galaxy is 
an indicator of an earlier formation time or equivalently an advanced stage of evolution for 
the group, greater excess H I content is expected, provided the diffuse intragroup gas can 
survive long enough compared to the mass assembly time. We find the mean and median 
for the excess gas mass fraction is 32 % and 23 % for the spiral dominated HCGs and the 
mean and the median for non-spiral dominated HCGs are 42 % and 49 %, respectively (see 
Table [7]). The graphical comparison shown on the bottom panel in Figure [12] shows the 
distribution of the spiral and non-spiral dominated groups. Most of the spiral dominated 
groups have lower excess gas mass fractions where as the non-spiral dominated groups cover 
the entire range. However, there is a some overlap between the two sub-samples. Based 
on the Analysis of Variance tesfl one can conclude that these two sub-samples are drawn 
from different distributions with more than 75% confidence. The weakness of the distinction 
between the two classes may indica te that the IG M evolution timescale is shorter than the 
mass assembly time ( ~ 4 Gyrs see lBarneslll989l ). 

Ano ther interesting observation base d on the VLA imaging of the non-spiral dominated 
systems (iBorthakur et al. in prep II2010I ) is that most of the H I mass is contributed by a 
single H I rich spiral galaxy, except for in the case of HCG 97 and HCG 37, where no H I 
was observed in the VLA imaging. In all the cases, the H I rich spiral galaxy is relatively 
undisturbed and shows no sign of tidal interaction. Interactions in galaxies can be seen 
in term s of disturbed H I morphology even when they are well away from the center of the 
cluster (jChung et al.ll2007l ). The most plausible explanation for a galaxy with an undisturbed 
H I distribution in a compact group is that it is new to the group and has not yet passed 
through the group center. So for all these systems, the VLA detected H I includes both the 
tidally stripped HSB gas as well as the H I mass of the undisturbed galaxy, thus diluting 
the excess gas fraction's ability to trace evolution of tidal debris. This may have caused the 
excess gas mass fraction in the non-spirals to span a much wider range as seen in Figure [T2J 



5. CONCLUSION 

A total of 26 HCGs have been observed in the 21 cm H I emission with the 100-meter 
GBT. We detected H I in almost all the compact groups with the exception of HCG 35. The 
derived H I masses ranged between 0.5 — 28.5 x 10 9 M . In HCG 40 and HCG 48, we 
detected 3 to 5 times the previously detected H I (compiled by VM2001) reducing their H I 
deficiencies significantly. 



5 An extension of Students t-test to more than two samples. 
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For our distance limited complete sample of 22 HCGs, we compared our GBT spectra 
with the corresponding VLA spectra for the region covered by the GBT beam. Strong 
differences between the two spectra can be seen for some groups while others showed very 
little or no differences at all. The total flux recovered by the GBT was always equal to 
or larger than that by the VLA (with the exception of HCG 40) as expected from correct 
calibration. Although the GBT had larger line flux, the shape of the two H I profiles for most 
of the groups are quite similar. There are fine differences in the form of broad faint emission 
and wing-like emission features that were only detected by the GBT. We have taken extra 
steps to ensure that these differences are not a result of calibration errors. Perfect match 
between the VLA and the GBT spectra of HCG 16, 93, and others supports our claim of 
calibration accuracy. By comparing the instrumental properties and data analysis techniques 
for the two instruments, we infer that this excess H I is in the form of a faint diffuse extended 
neutral IGM. In a few groups, spatial filtration of extended H I structures of angular size 
>12-15' by the VLA is also plausible. 

Owing to their high spatial density and low velocity dispersion, tidal interactions are 
ubiquitous in HCGs. This results in a substantial deposition of cold gas from the ISM of 
individual galaxies to the hot IGM. The detection of neutral IGM in HCGs confirms that the 
cold tidally stripped filaments and clouds survive the conductive heating and expansion for 
longer than a typical group dynamical timescale (> 400 Myrs). Based on our calculations, 
we constrain the radius of an H I cloud to be > 200 pc to survive over the dynamical time of 
the group. Furthermore, the presence of neutral medium also constrains the temperature of 
the IGM to be lower than 4 x 10 6 K. We also conclude that the diffuse neutral gas content 
should increase with interactions between member galaxies and hence is related to the age 
of the group. 

We characterize the excess gas content in terms of the total H I content in the galaxy- 
group. The excess gas mass fraction, /= M excess /M GBT , varies from almost none to 81%. 
We discovered a much larger fraction of excess gas (median=51 %) in slightly H I deficient 
groups as compared to H I normal groups ( median=26 %) or highly H I deficient groups 
(median=20 %). For groups with H I deficiency less than 0.4 (= 2 a of deficiency estimation), 
the excess gas fraction increases with increasing deviation from H I normalcy. The opposite 
trend is seen in the case of high-surface brightness H I detected by the VLA. We conclude 
that the diffuse neutral gas represents an intermediate stage in the evolution of tidal debris 
in the hot IGM between HSB H I and ionized hydrogen. 

We also examined the excess gas content in terms of group evolution based on the 
scenario proposed by VM2001. According to this evolutionary scenario, the tidal interactions 
increasingly deposit ISM from individual galaxies into the IGM as a group evolves and hence 
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based on the H I morphology one can predict the stage of evolution. The tidal debris, which 
survives the hot IGM, will evolve into the diffuse neutral IGM and is expected to be a 
substantial fraction of the total H I content. Thus the excess H I content is expected to be 
correlated with the evolutionary phase. From our analysis we found that this is indeed the 
case. The groups in the final phase of evolution had a larger amount of excess H I (average 
of 48%) as compared to groups in phase 1 (21%) and phase 2 (30%). Furthermore, we find 
a clustering in the distribution with only two systems showing the transition phase. This 
suggests that the transition is rapid and hints at the presence of a physical threshold like a 
characteristic age or an interaction threshold aiding the transition. We also examined the 
dependence of excess gas content on the Hubble type of the central dominant galaxy in the 
group. While spiral dominated groups tend to be have a lower excess gas mass fraction, the 
non-spiral dominated groups have a wide range of values covering the entire range. 

The existence of cold neutral IGM in HCGs offers a strong support for the survival of the 
tidally stripped ISM from individual galaxies. A significant amount of H I is tidally stripped 
from the galaxies as they interact with the group potential and other members and ultimately 
ends up in the diffuse IGM. Thus, the excess gas content gives valuable information on the 
evolutionary history of groups. In addition, the neutral gas distribution also holds the key to 
understanding current and future star formation in groups, which in-turn can considerably 
change the morphology and spectral properties of the galaxies. Cooler regions of the IGM 
might form clumps with sufficient densities to fall back into the galaxies, and providing them 
with fresh fuel for star-formation. The presence of neutral IGM also raises the possibility of 
non-disk star formation. Therefore, this discovery opens up a host of unanswered questions 
and we hope that a GBT-VLA combined map will provide clues to some of these questions. 
Sensitive X-ray masses and maps along with a UV absorption study of the IGM would be 
the next crucial step in bringing out the full picture. 

This manuscript has benefited from insightful discussions and comments by J. Hibbard, 
J. Kwan, F. Lockman, A. Minter, H. Mo, T. Ponman, J. Rasmussen, T. Tripp, J. van 
Gorkom, and K. C. Xu. The authors are grateful to the observatory staff at the GBT 
who made these observations possible. Support for this work was provided by the NSF 
through award GSSP06-0002 from the NRAO. LVM is partially supported by DGI Grant 
AYA 2005-07516-C02-01 and Junta de Andaluca (Spain). This research has made use of the 
NASA/IPAC Extragalactic Database (NED) operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 
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Fig. 1. — H I profile of the 26 Hickson Compact Groups observed with the 100 meter 
diameter GBT. The redshifts of individual members are marked with short vertical lines on 
the spectra. 
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Fig. 1. — continued. 




Fig. 1. — continued. 
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Observed HI Mass (Log M HI ) 

Fig. 2. — The H I mass observed using the GBT is compared with the predicted H I mass 
for each of the 22 HCGs from our complete sample as well as the 4 additional sources. The 
diagonal line drawn indicates H I normalcy. A clear majority of the groups appear on the 
H I deficient side. Typical uncertainties in both masses are shown using error bars on the 
bottom right corner. The additional sources are H I normal. No H I was detected in HCG 35 
and consequently an uppper limit is presented. 
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Deficiency 



Fig. 3. — Histogram showing Hi-deficiency evaluated using our GBT measurements. The 
dashed line represents zero deficiency and the dotted line represents the expected distribution 
for an H I normal sample (see § I3.2p . Typical uncertainty in Hi-deficiency is estimated to 
be about 0.2 in the log scale. HCG 35 is not included in this plot since no H I was detected 
and hence reliable H I deficiency value could not be estimated. 
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Fig. 4. — Comparison of the GBT spectrum (gray line) with VLA spectrum (dark line) for 
our complete sample of 22 groups as well as two of the addition sources- HCG 26 and 48 
in the last two panels. The VLA spectra presented here were obtained by superimposing 
the GBT beam on the VLA imaere (see 8 13.31). Details of the VLA data, are summarized in 
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Fig. 4. — continued. 
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Fig. 4. — continued. 
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Fig. 4. — continued. 
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Log(M HI ) Log(M HI ) 

Fig. 5. — Histogram showing the distribution of total HI mass detected by the GBT (left) 
and excess HI mass (right) for our complete sample with the exception of HCG 40. On an 
average the excess gas is more than a third of the total H I detected by the GBT, although 
individual values vary from 5% to 81% within the sample. 
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Fig. 6. — An overlay of the VLA spectrum (solid black) and the GBT spectrum in (black 
dash) obtained by binning GBT data to VLA spectral resolution for HCG 92. The high 
velocity wings associated with the peak at 6620 km s _1 extends from 6200 to 6550 km s _1 . 
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Fig. 7. — An overlay of the VLA spectrum in solid black obtained using the Multi-scale 
CLEAN (MSC) algorithm in AIPS over the GBT and the VLA spectra for HCG 100. The 
VLA MSC spectrum recovers 20% more flux than the VLA spectrum although it does not 
fully recover the GBT flux. This confirms that there is flux at shorter baselines or larger 
scales in the VLA data, although the VLA misses the remainder of the GBT flux due to the 
absence of shorter baselines (or zero baseline) in the data. This translates to a spatial scale 
of ~15' for the H I structures on the plane of the sky. 
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Fig. 9. — Plots showing cloud evaporation timescale regimes as a function of cloud radius 
and density of the IGM for different sets of IGM temperature and cloud density. The Tigm 
and n c i values are printed on the top of each plot. All notations are same as Figure 
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Fig. 10. — The diffuse gas fraction is plotted as a function of the GBT beam size (in terms of 
the group radius). On the right is a histogram HCGs is shown as a function of the diffuse gas 
fraction. The plots do not include the extremely compact group HCG 79 (GBT beam/group 
diameter ~7) with 21% of excess gas. 
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Fig. 11. — The diffuse gas fraction is plotted as a function of the GBT beam size in kilopar- 
secs. 
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Fig. 12. — Comparisons of the diffuse gas mass fraction for our complete sample (except 
HCG 40) in terms of H I deficiency (top), group evolutionary phase (middle), and the 
dominant group galaxy type (bottom). The bimodal distribution of the diffuse gas mass 
fraction is clearly seen in these histograms. 
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Fig. 13. — Variation of diffuse gas content and the high surface brightness (HSB) VLA 
detected H I in terms of the predicted H I mass for 11 HCGs with H I deficiency < 0.4. The 
dotted line marks the case of H I normalcy. 



Table 1. Information on the GBT Sample. 
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5.0 
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c Evaluated using same cosmology as VM2001 - Q m = 1 and Q\ = 



jde Carvalho et al.l l ll997l) 
IVerdes-Monteneero et al.l d2005h 
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Table 2. Results from our GBT observations and comparison with previously published data. 



Source 


HI mass 
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0.20 
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0.28 


0.50 
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0.5 
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9.89 


10.03 


-0.14 


-0.14 
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10.23 


10.45 


-0.22 


-0.19 
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HCG 48 


1.2 


0.5 
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Hfficksonl dl982h 

b Evaluated using same cosmology as VM2001 - f2 m = 1 and Q\ = 

c Noise measured for smoothed data of resolution 10 km s _1 H I mass estimated assuming a line width A V=100 km s 



d Within the GBT beam. Values have been borrowed from VM2001. 

e Calculated by the authors using the same procedure as VM2001 taking into account the GBT beam (~ 9.1') coverage. 

f GBT beam is smaller than the extent of H I in VLA maps and hence values presented are significantly lower than total H I content. 

g No H I was detected. Values refer to mass corresponding to 3cr noise level assuming a line-width of 600 km s _1 . 



Table 3. A summary of our VLA data (used for comparison with GBT spectra). 



Source 


Array 


Synth. Beam 


Channel Width 
km s — 1 


Line width 


Reference 


HCG 7 


D 
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21.2 
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VLA archive 


HCG 10 
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x 53" 


21.3 
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Borthakur et al. in prep (2010) 
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1) 
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1) 
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21.5 




Yun et al. m prep (2010) 
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21.5 


550 


Yun et al. in prep (2010) 
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21.4 
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Williams et al. (2002) 
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60' 
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21.4 


600 


Borthakur et al. in prep (2010) 


HCG 97 


D 


69' 


x 55" 


21.6 
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Borthakur et al. in prep (2010) 


HCG 100 


D 


60' 


x 56" 


21.4 


300 


Borthakur et al. in prep (2010) 


HCG 26 


C 


26' 


x 18" 


11.0 


500 


Williams & Van Gorkom (1995) 


HCG 48 


CD 


53" 


x 39" 


21.0 


150 


Yun et al. in prep (2010) 



a Only part of the single dish emission was covered by the VLA spectral range 
b Sum of multiple emission regions 



Table 4. Filling factor calculation for the excess gas seen in the GBT spectra. 
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HCG 67 


15 


6.9 


I 


0.53 


42 


1.50 


11 


HCG 68 


54 


30.3 


1 


0.31 


42 


1.24 


566 d 


HCG 79 


21 


9.0 


II, III 


0.90 


42 


19.64 


3 


HCG 88 


21 


24.6 


I 


0.31 


42 


9.38 


10 


HCG 90 


44 


2.3 


III 


0.34 


12 


1.41 


32 


HCG 91 


38 


86.9 


1 


0.51 


12 


0.89 


254 d 


HCG 92 


51 


86.3 


I, II, III 


0.21 


42 


7.09 


39 


HCG 93 


8 


2.2 


I 


0.33 


42 


1.20 


9 


HCG 97 


55 


23.5 


I, III 


0.52 


12 


1.60 


45 


HCG 100 


55 


49.3 


I, II 


0.30 


42 


1.10 


209 d 


HCG 26 


28 


79.9 


1 


0.54 


42 


14.70 


8 


HCG 48 


11 


5.1 


I 


0.42 


42 


2.50 


34 



a Excess H I fraction, M excess 

/M G BT 

b See §rO]for details. 

c Spectrum used for comparison is incomplete in its velocity coverage. 

Significant spatial filtration of H I by the VLA in groups with the H I structures are of the order 10-15'. 



Note. — HCG 40 is not included in the above table as excess gas mass could not be determined precisely due to uncertainty in calibration 
of the GBT data for this group. 



Table 5: The excess gas mass fraction for the three H I deficiency classes defined in § 13.21 



Normal H I content 


Slightly H I deficient 


Highly H I deficient 


Source 


% detection 


Source 


% detection 


Source 


% detection 




^-excess f^-GBT 




M excess /Mgbt 




Mezcess/Mc-BT 


HCG 25 


15 % 


HCG 16 


5 % 


HCG 93 


8 % 


HCG 79 


21 % 


HCG 58 


13 % 


HCG 67 


15 % 


HCG 23 a 


30 % 


HCG 10 


17% 


HCG 7 


16 % 


HCG 68 


54 % 


HCG 88 


21 % 


HCG 30 


23 % 






HCG 91 


38 % 


HCG 90 


44 % 






HCG 92 


51 % 


HCG 44 


60 % 






HCG 97 


55 % 










HCG 100 


55 % 










HCG 15 


56 % 










HCG 31 


70 % 










HCG 37 


81 % 






Average 


30 % 


Average 


42 % 


Average 


28 % 


Median 


26 % 


Median 


51 % 


Median 


20 % 



"Spectrum used for comparison is incomplete in its velocity coverage. 
Note. — HCG 40 is not included in the analysis (see § 14.41 for details). 



-56 - 



Table 6: The excess gas mass fraction for the three phases of group evolution 



Phase 1 




Phase 2 


Phase 3 


Source 


% detection 


Source 


% detection 


Source 


% detection 




M e xcess 


/M GBT 




M excess /McBT 




^-excess /^-GBT 


HCG 67 


15 % 


HCG 16 


5% 


HCG 93 


8 % 


HCG 7 


16 


% 


HCG 58 


13 % 


HCG 30 


23 % 


HCG 88 


21 


% 


HCG 25 


15 % 


HCG 90 


44 % 


HCG 79 


21 


% 


HCG 10 


17% 


HCG 92 


51 % 


HCG 23 b 


30 


% 


HCG 91 


38 % 


HCG 68 


54 % 








HCG 100 


55 % 


HCG 97 


55 % 








HCG 31 


70 % 


HCG 15 


56 % 












HCG 44 


60 % 












HCG 37 


81 % 


Average 


21 


% 


Average 


30 % 


Average 


48 % 


Median 


21 


% 


Median 


17% 


Median 


54 % 



a VM2001 

& Spectrum used for comparison is incomplete in its velocity coverage. 

Note. — HCG 40 is not included in the analysis (see § 14.41 for details). 
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Table 7: The excess gas mass fraction for the spiral dominated and non-spiral dominated 
HCGss a . 

S : )i >';■ I Non-spiral 



Source % detection Source % detection 

M excess /M GBT M 

excess 

/M GBT 



HCG 16 


5 


% 


HCG 93 


8 ( 


Vo 


HCG 58 


13 


% 


HCG 67 


15 


% 


HCG 25 


15 


% 


HCG 79 


21 


% 


HCG 7 


16 


% 


HCG 90 


44 


% 


HCG 10 


17 


% 


HCG 68 


54 


% 


HCG 88 


21 


% 


HCG 97 


55 


% 


HCG 30 


23 


% 


HCG 15 


56 


% 


HCG 23 b 


30 


% 


HCG 37 


81 


% 


HCG 91 


38 


% 








HCG 92 


51 


% 








HCG 100 


55 


% 








HCG 44 


60 


% 








HCG 31 


70 


% 








Average 


32 % 


Average 


42 


% 


Median 


23% 


Median 


49 


% 



iHicksonl (119821 ) 

& Spectrum used for comparison is incomplete in its velocity coverage. 



Note. — HCG 40 is not included in the analysis (see § 14.41 for details). 



